The fate of carbon-based life as a function of the light quark mass 
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The Hoyle state plays a crucial role in the helium burning of stars that have reached the red giant stage. The 
close proximity of this state to the triple-alpha threshold is needed for the production of carbon, oxygen, and 
other elements necessary for life. We investigate whether this life-essential condition is robust or delicately fine- 
tuned by measuring its dependence on the fundamental constants of nature, specifically the light quark mass and 
the strength of the electromagnetic interaction. We show that there exist strong correlations between the alpha- 
particle binding energy and the various energies relevant to the triple-alpha process. We derive limits on the 
variation of these fundamental parameters from the requirement that sufficient amounts of carbon and oxygen 
be generated in stars. We also discuss the implications of our results for an anthropic view of the Universe. 
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Life as we know it depends on the availability of carbon and 
oxygen. These two essential elements are produced during he- 
lium burning in red giant stars. The initial reaction is the so- 
called triple-alpha process, where three helium nuclei fuse to 
generate 12 C. This can be viewed as a two-step process. First, 
two 4 He nuclei combine to form an unstable, but long-lived 
8 Be resonance. This 8 Be resonance must then combine with a 
third alpha-particle to generate carbon. By itself, this process 
cannot explain the observed abundance of carbon in the Uni- 
verse. Therefore, Hoyle postulated that a new excited state 
of 12 C, a spinless even-parity resonance near the 8 Be-alpha 
threshold, enhances the reaction JlJ. Soon after this predic- 
tion, the new state was found at Caltech J2, 0] and has since 
been investigated in laboratories worldwide. The measured 
energy of this second + state is e = 397.47(18) keV above 
the triple-alpha threshold, while the total and radiative widths 
are known to be T tot = 8.3(1.0) eV and T y = 3.7(5) meV, 
respectively. The reaction rate for the (resonant) triple- alpha 
process is approximately given by |01 



r 3a « T 7 {Njk B Tf exp(- £ /fc B T), 



(1) 



with N a the alpha-particle number density, T the stellar tem- 
perature and k B Boltzmann's constant. Due to the exponen- 
tial dependence, e is the dominant control parameter of this 
reaction. Here, we study the dependence of e upon the fun- 
damental parameters of the strong and electromagnetic (EM) 
interactions. 

Given its role in the formation of life-essential elements, the 
Hoyle state has been called the "level of life" Jfl (see Ref. |0] 
for a thorough discussion of the history of this issue). Thus, 
it is often considered a prime example of the anthropic princi- 
ple, which states that the observed values of the fundamental 
physical and cosmological parameters are restricted by the re- 
quirement that life can form to observe them, and that the cur- 
rent Universe be old enough for that to happen 10, HI ■ In the 



context of cosmology and string theory, consequences derived 
from anthropic considerations have had considerable impact 
(see e.g. Refs. HEl). 

Several numerical studies have investigated the impact of 
changes in the Hoyle state energy. Livio et al. Ill 111 modified 
the value of e by hand and performed calculations involving 
the triple-alpha process in the core and helium shell burning 
of helium up to the asymptotic giant branch stage in stellar 
evolution. They concluded that a ~ 60 keV change in e could 
be tolerated, and thus the amount of fine-tuning required was 
not as severe as first believed. 



A more microscopic calculation was performed by Ober- 
hummer et al. iR FlHi in terms of a nuclear cluster model based 
on a simple two-nucleon (NN) + EM interaction. This NN 
interaction was formulated in terms of one strength parame- 
ter, adjusted to give a fair description of a-a scattering and 
the spectrum of 12 C. By modifying this coupling strength and 
the EM fine structure constant a em , the effect on carbon and 
oxygen production was analyzed. Outside of a narrow win- 
dow of ~ 0.5% around the observed strong force and ~ 4% 
around the observed Coulomb force, the stellar production of 
carbon and/or oxygen was found to be reduced by several or- 
ders of magnitude. However, this model of the strong force is 
not readily connected to the fundamental theory of the strong 
interactions, quantum chromodynamics (QCD), and its fun- 
damental parameters, the light quark masses. Therefore, it is 
not obvious how to translate the findings of Ref. 1U2I1 into an- 
thropic constraints on fundamental parameters. In this study, 
we shall address this pertinent question: What changes in the 
quark masses and the EM fine structure constant are consistent 
with the formation of carbon-based life? 

Over the last few years, we have developed a new method 
to study atomic nuclei and their properties from first princi - 
pies, termed nuclear lattice simulations. The key ingredients 
in this approach are, on the one hand, the chiral effective field 
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theory (EFT) of nuclear forces and, on the other hand, large- 
scale lattice Monte Carlo methods. The latter are also fruit- 
fully used in many other fields of science. Chiral nuclear EFT 
was introduced by Weinberg IU3I1 (for a first numerical imple- 
mentation, see Il4ll ) as a systematic tool to explore the conse- 
quences of spontaneous and explicit chiral symmetry breaking 
of QCD in a rigorous manner. The basic degrees of freedom 
are pions and nucleons, where the pions and their interactions 
carry the basic information of the chiral symmetry properties 
of QCD. In particular, one finds M% ~ (m u + m d ), so that 
any dependence on the light quark masses m u and m d can 
be translated into a corresponding dependence on the pion 
mass M w . In what follows, only the average light quark mass 
m q = (m u + rn d )/2 will be considered, as the effects of 
strong isospin violation due to m u 7^ m d are greatly sup- 
pressed for the reactions considered here. Chiral nuclear EFT 
is based on an order-by-order expansion of the nuclear po- 
tential. In this scheme, two-, three- and four-nucleon forces 
arise naturally, and their observed hierarchy is also explained. 
The nuclear forces have been worked out to high precision 
and applied successfully in few-nucleon systems for binding 
energies, structure, and reactions. For recent reviews, see 
Refs. lfl5i [l6ll . Within chiral nuclear EFT, the quark mass 
dependence of light nuclei and its impact on big ban g nu cle- 
osynthesis has already been studied, see e.g. Refs. 1 17142211 
and M23I1 for a related study. 

Monte Carlo simulations have been used to solve the nu- 
clear A-body problem (with A the atomic number) based on 
a lattice formulation [24-]. The lattice spacing of this dis- 
cretized space-time serves as an ultraviolet regulator. The nu- 
cleons are placed on the lattice sites, and the interactions are 
represented by pionic and (suitably chosen) auxiliary fields. 
Our periodic cubic lattice has a spacing of a = 1.97 fm and 
a length of L = 11.82 fm. In the time direction, our lat- 
tice spacing is a t = 1.32 fm, and the propagation time L t 
is varied in order to extrapolate to L t — > 00. The energies 
of the ground and excited states are obtained using projection 
Monte Carlo techniques ll25il26ll . More precisely, we compute 
Z A (t) = (ip A \ cxp(—Ht)\ip A ) for a given A-nucleon system 
at large Euclidean time t in order to extract the energies of the 
low-lying states, see also Ref. 12711 for more details. 

The leading order (LO) contribution to the NN force 
emerges from the one-pion exchange potential (OPEP) and 
(smeared) 5-wave contact interactions. This improved LO 
action forms the basis of our projection Monte Carlo simu- 
lations, while all higher-order terms including the Coulomb 
interaction, corrections to the NN force and three-nucleon 
forces, are treated in perturbation theory. All parameters of 
H are fixed from two- and three-nucleon data, enabling pre- 
dictions for all heavier nuclei. So far, such calculations have 
been performed up to next-to-next-to-leading order (NNLO), 
achieving a good description of nuclei up to A = 12. 

We have performed the first ab initio calculations for the 
energy |25|] and structure of the Hoyle state 112 611 using the 
nuclear lattice formalism (see Ref. 112 811 for a no-core shell 
model calculation of the spectrum of 12 C employing chiral 



EFT forces). In our approach, the hadronic interactions of 
the nucleons with themselves and with pions can be modified 
easily. Our analysis of the dependence upon the strength of 
the Coulomb interaction is therefore straightforward. For the 
dependence on m q , we also need information about the quark 
mass dependence of the hadronic interactions. In turn, such 
dependences can be given as a function of M„. 

We shall restrict ourselves to values of near the phys- 
ical point, with {SM^/M^l < 10%. Such small changes can 
be treated in perturbation theory. The M^-dependence of the 
OPEP and the nucleon mass m N is determined in chiral per- 
turbation theory utilizing constraints from lattice QCD, see 
Ref. I3H1 for more details. To retain model independence, we 
do not rely on the chiral expansion of the NN contact interac- 
tions. Instead, we express our results in terms of the deriva- 
tives of the inverse spin-singlet and spin-triplet NN scattering 
lengths with respect to the pion mass, 



A. 



-4, 



(2) 



which parameterize the A/„.-dependence of the short-range 
nuclear force and can be measured in lattice QCD. We do not 
consider Af^-dependent short-range effects beyond the ones 
introduced above. 1 Thus, the variation of a given nuclear en- 
ergy level E i takes the form 



dE, 



dE, 



dAL 



with Xi 



dM OPE 
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(3) 
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3E l 
dC n 



dE, 



cf 



dm N /dMJ mh , x 2 = dg^ N /dM n \ M ^ etc., 
where X ph denotes the value of X for the physical M v . The 
terms in Eq. © represent different contributions to the pion 
mass variation. First, there is the explicit dependence on M T 
through the pion propagator in the OPEP. Second, we include 
the dependences on M w through the nucleon mass m N and 
g N = g A j (2F T ), with g A the nucleon axial-vector coupling 
and F„ the weak pion decay constant. Finally, we have the 
AT^-dependences from the strengths of the NN contact inter- 
actions C and Cj, which are expressed through the deriva- 
tives given in Eq.©. Therefore, the problem reduces to the 
calculation of various derivatives of the nuclear energy levels 
using lattice Monte Carlo techniques and the determination of 
the coefficients x 1 ...x±. The derivatives of E i in Eq. <[3j are 
computed by evaluating the expectation value of the deriva- 
tive of the lattice Hamiltonian H with respect to A/° PE , m N , 
g nN , C and Cj. This involved the generation of O(10 7 ) sta- 
tistically independent pion- and auxiliary field configurations 



1 The correlations observed for various energy differences (as discussed be- 
low) indicate that the dynamics of interest are largely governed by the large 
5-wave NN scattering lengths. Higher-order M n -dependent short-range 
terms are therefore expected to play a minor role. 
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FIG. 1 : Correlation of the a-particle energy E 4 with the energy dif- 
ferences pertinent to the triple-alpha process under variation of 4>,t 
in the range { — 1 ... 1}. 



on the Blue Gene/Q supercomputer JUQUEEN using the hy- 
brid Monte Carlo algorithm. The explicit form of H can be 
found in Ref. JH. 

The values of x 1 and x 2 can be obtained from lattice QCD 
combined with chiral extrapolations (see e.g. Ref. l3(ill for a 
recent review on lattice QCD and determinations of the nu- 
cleon mass variation). We exchange x 3 and x 4 for A s and 
A t by consideration of the A/^-dependence of NN scattering 
in a cubic box. We may then compute the energy differences 



AE h 



ET 



E a 



E 4 and AE b 



E 8 - 2E 4 , where ET 2 



is the energy of the Hoyle state and E 4 8 the ground-state en- 
ergies of the 4 He and 8 Be nuclei, respectively. Note also that 
We find 



e ee AE h + AE b 



dAE h 



dAE h 



de 



-0.455(35)4 - 0.744(24)4 + 0.051(19), 



-0.117(34)4, - 0.189(24)A t + 0.013(12), 



dM„ 



= -0.572(19)4 - 0.933(15)4 + 0.064(16), 



(4) 



where the parentheses represent the one-standard-deviation 
stochastic (Monte Carlo) and extrapolation error, combined 
with the uncertainty in x 4 2 (as explained in Ref. If3l"tl ) which 
affects only the constant (OPEP) terms above. 

The results in Eq. are intriguing. First, we note that 
(dAEjdM„)/(dAE b /dM % ) ~ 4, thus AE h and AE b can- 
not be independently fine-tuned. Such behavior can readily 
be explained in terms of the a-cluster structure of the Hoyle 
state and 8 Be. Further correlations are visualized in Fig. [T] 
where the relative changes in AE bl AE h and e are shown as 
a function of relative changes in the ground state energy E 4 
of the a-particle. We define K\ = (dX/dM w )M n /X as the 
relative variation of X with respect to M n . Fig. [T] provides 



clear evidence that the alpha binding energy is strongly cor- 
related with AE h , AE h , and e. Such correlations related to 



carbon production have been speculated upon earlier Hi UI32H . 
Second, we note that there is a special value for the ratio of 
4 t° 4, given by 



-1.5 



(5) 



where the pion mass dependence of AE h , AE b , and e be- 
comes small (compared to the error bars). 

We have expressed all our results in terms of the quantities 
4 t> the quark mass dependence of which was considered at 



next-to-leading order (NLO) in Ref. 11811 . with a recent up- 
date to NNLO HI. That analysis gives A s = 0.29toJ;j and 
where the errors reflect the theoretical un- 



A t = 



-0.18 



+0.10 
-0.10' 



certainties. As expected, these values of A s t are of natural 
size. Taking into account correlations in the calculation of 
4t> we f 111 ^ A s /A t = —1.6^{®. Interestingly, the central 
value is very close to the result given in Eq. (|5), for which the 
pion mass dependences of AE h , AE b , and e are all approx- 
imately zero (within error bars). In the future, a reduction of 
the uncertainty in A s t is desirable. This can be addressed by 
lattice QCD calculations of NN systems. For recent studies, 
seeRefs. SB- 

We now use the reaction rate in Eq. (Q3 to draw conclu- 
sions about the allowed variations of the fundamental con- 



stants. From the stellar modeling calculations in Ref. Ill 211 . we 
find that sufficient abundances of both carbon and oxygen can 
be maintained within an envelope of ±100 keV around the ob- 
served value of e. Allowing for a maximum shift of ±100 keV 
in e translates into bounds on the variations of m q . In Fig. [2] 
we show "survivability bands" for carbon-oxygen based life 
due to 1% and 5% changes in m q (in terms of A s and A t ). 
To be precise, for a 5% change in m„, A t must assume values 
within the red (narrow) band to allow for sufficient produc- 
tion of carbon and oxygen. The most up-to-date knowledge 
of these parameters is depicted by the data point with hori- 
zontal and vertical error bars. This NNLO determination of 
4 t shows that carbon-based life survives at least a ~ 0.7% 
shift in m q . In addition to this "worst-case scenario", we find 
that the theoretical uncertainty in A s t is also compatible with 
a vanishing de/dM„ (complete lack of fine-tuning). Given 
the central values of A st , we conclude that variations of the 
light quark masses of 2 — 3% are unlikely to be catastrophic 
to the formation of life-essential carbon and oxygen. 

We may also compute the corresponding changes induced 
by variations of the EM fine-structure constant a om . On the 
lattice, the EM shift receives contributions from the long- 
range Coulomb force and a short-range proton-proton con- 
tact interaction. The latter contains an unknown coupling 
strength, which allows for the regularization of QED on the 
lattice. We have fixed its finite part from the known EM con- 
tribution to the a-particle binding energy. The dependence 
of the E i on a em can then be calculated. By expressing 
the EM shifts as (dX/da cm )\ a ^ ~ Q(X)/a em , we find 
Q{AE b ) = 1.19(8) MeV, Q(AE h ) = 2.80(10) MeV and 
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FIG. 2: "Survivability bands" for carbon-oxygen based life from 
Eq. I0, due to 1% and 5% changes in m, in terms of the param- 
eters A B and A t . The most up-to-date NNLO analysis of A Stt is 
depicted by the data point with horizontal and vertical error bars. 

Q(s) = 3.99(9) MeV. For fixed m„, a variation of a cm by 
±100 keV/Q(e) « 2.5% would thus be compatible with the 
formation of carbon and oxygen in our Universe. This is con- 
sistent with the ~ 4% bound reported in Ref. fl. 

In summary, we have presented ab initio lattice calculations 
of the dependence of the triple-alpha process upon the light 
quark masses and the EM fine structure constant. The posi- 
tion of the 8 Be ground state relative to the 2-a threshold, as 
well as that of the Hoyle state relative to the 3-a threshold, 
appears strongly correlated with the binding energy of the a- 
particle. We also find that the formation of carbon and oxy- 
gen in our Universe would survive a change of ~ 2% in m„ 
or ~ 2% in a em . Beyond such relatively small changes, the 
anthropic principle appears necessary at this time to explain 
the observed reaction rate of the triple-alpha process. In order 
to make more definitive statements about carbon and oxygen 
production for larger changes in the fundamental parameters, 
a more precise determination of A s and A t is needed from 
future lattice QCD simulations. 
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